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We present measurements of the energy dependence of event-by-event fluctuations in the K/tt 
and (p + p)/7r multiplicity ratios in heavy ion collisions at the CERN SPS. The particle ratio 
fluctuations were obtained for central Pb+Pb collisions at five collision energies, ^/s NN , between 6.3 
and 17.3 GeV. After accounting for the effects of finite-number statistics and detector resolution, 
we extract the strength of non-statistical fluctuations at each energy. For the K/tt ratio, larger 
fluctuations than expected for independent particle production are found at all collision energies. 
The fluctuations in the (p + p) /tt ratio are smaller than expectations from independent particle 
production, indicating correlated pion and proton production from resonance decays. For both 
ratios, the deviation from purely statistical fluctuations shows an increase towards lower collision 
energies. The results are compared to transport model calculations, which fail to describe the energy 
dependence of the K/tt ratio fluctuations. 



I. INTRODUCTION 

Quantum Chromodynamics predicts that at suffi- 
ciently high temperature, strongly interacting matter will 
undergo a phase transition from hadronic matter to a 
state characterized by quark and gluon degrees of free- 



dom, the quark-gluon plasma (QGP) [l|. Experimen- 
tally, strongly interacting matter under extreme condi- 
tions can be created in heavy ion collisions at highly rela- 
tivistic energies. Experiments have been performed over 
a very large range of center of mass collision energies, 
y/s NN7 including 2.3 GeV < y/s NN < 4.9 GeV at the 



Brookhavcn Alternating Gradient Synchrotron (AGS), 
6.3 GeV < y/s^ < 17.3 GeV at the CERN Super Proton 
Synchrotron (SPS) and 19.6 GeV < ^/s^ < 200 GeV at 
the Brookhaven Relativistic Heavy Ion Collider (RHIC). 

In the SPS energy range, several hadronic observables 
in central Pb+Pb collisions show qualitative changes in 
their energy dependence. The ratio of average (K + ) to 
(7r + ) yields exhibits non-monotonic behavior in the low- 
energy SPS range, close to y/s NN w 7.6 GeV 0. In 
the same energy range, the slopes of the hadron trans- 
verse momentum distributions show an approximately 
constant value, after a rapid rise at lower energies [3j. 
These features are not observed in elementary interac- 
tions and appear to be unique characteristics of heavy ion 
collisions. Moreover, these observations have not been 
reproduced in hadronic transport models or statistical 
models of hadron production. The data are consistent 
with the expected signals of the onset of a phase tran- 
sition in heavy ion collisions at low SPS energies 0, E|- 
Recent measurements in Au+Au collisions at the highest 
available energies at RHIC indicate that in these colli- 
sions a highly collective, nearly thermalized system is 
formed although the nature of the underlying de- 
grees of freedom is still under debate. A comparison of 
SPS to RHIC data shows a rather smooth evolution of 
hadronic observables from the higher SPS energies up to 
the highest RHIC energies d, @ . 

Further information about the existence and nature 
of a phase transition in the SPS energy range can pos- 
sibly be gained from studies of event-by-event fluctua- 
tions of final state hadron distributions. Several mecha- 
nisms have been proposed that could lead to such event- 
by-event fluctuations, including overheating-supercooling 
fluctuations due to a first order phase transition with 
large latent heat 0, and fluctuations due to a phase of 
coexisting confined and deconfined matter (mixed phase) 
in varying relative abundances [f|. The presence or ab- 
sence of such fluctuations might contain otherwise in- 
accessible information about the nature and order of a 
phase transition at these energies Q . 

Additional motivation for studies of fluctuations in 
kaon and proton production comes from theoretical stud- 
ies of the QCD phase diagram d, [l(| [H[. ^ has been 
argued that the transition from hadronic to deconfined 
matter changes from a first order transition at large 
baryo-chemical potential, (j,b, to a second-order transi- 
tion or a cross-over at small [Lb- This implies the exis- 
tence of a critical endpoint at the end of the first order 
transition line. If freeze-out of a bulk strongly interact- 
ing system occurred close to the critical point, observable 
non-statistical fluctuations in particle momentum distri- 
butions or baryon number correlations could be induced. 
While theoretical predictions for the location of the criti- 
cal point and the strength of expected fluctuation signals 
are still under debate, it is clear that a positive experi- 
mental signal for fluctuations related to the critical point 
could lead to major progress in our understanding of the 
QCD phase diagram. 



To obtain experimental information directly related to 
the nature of the phase transition from fluctuations, sys- 
tematic measurements in various regions of the /j b versus 
temperature T phase diagram are necessary. The en- 
ergy density, temperature and baryo-chemical potential 
probed in the reaction can be controlled by varying the 
incident energy of the colliding system. This will allow us 
to test if the observed structure in the energy dependence 
of event-averaged hadronic signatures is reflected in the 
collective phenomena probed by event-by-event fluctua- 
tions. Preliminary results by NA49 showing a significant 
energy dependence of fluctuations in the kaon to pion 
and proton to pion multiplicity ratios were shown in [l2| ■ 
In this paper, we present final results on the energy de- 
pendence of event-by-event fluctuations of these ratios in 
central Pb+Pb collisions over the SPS energy regime. 

Throughout the paper, we will define the event-by- 
event kaon to pion ratio as Kjix = jf £ ~fN K ~ and 
the event-by-event proton to pion ratio as (p + p)/tv = 
N '+ N P ■ Here, N^±, N K ±, N p and N p denote the 
observed pion, kaon, proton and anti-proton yields, re- 
spectively, within the acceptance for a given event. It 
is important to note that all multiplicities, particle ra- 
tios and fluctuations reported in this paper refer to the 
observed yields in the experimental acceptance region. 
For the fluctuation measurement, an extrapolation to full 
phase space would only be possible using an assumption 
about the (unknown) nature of the underlying correla- 
tions or fluctuations. To allow a detailed comparison of 
theoretical calculations to the results presented here, 3- 
dimensional acceptance tables for each collision energy 
are available online [l3[. 



II. EXPERIMENTAL SETUP AND DATA 
SELECTION 

The data presented here were taken with the NA49 ex- 
periment during runs from 1996 to 2002. The NA49 setup 
is shown schematically in Fig. [1] and described in detail 
in [l4j. Particle trajectories are measured using four large 
volume Time Projection Chambers (TPCs). Two TPCs, 
VTPC-1 and VTPC-2, are placed in the magnetic field of 
two super-conducting dipole magnets. Two other TPCs, 
MTPC-L and MTPC-R, are positioned downstream of 
the magnets. The MTPCs were optimized for a high pre- 
cision measurement of the specific ionization energy loss 
dE/dx which provides the main method of particle iden- 
tification for this analysis. The target, a thin lead foil of 
about 0.01 interaction length for Pb ions, was positioned 
80 cm upstream from VTPC-1. 



A. Event selection 

At all five beam energies the most central Pb+Pb col- 
lisions were selected based on the energy E veto deposited 
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FIG. 1: (Color online) The experimental set-up of the NA49 experiment [1J| and target arrangement (a) for Pb+Pb collisions. 

TABLE I: Details of the data sets used in this analysis. For all data sets, the 3.5% most central events were selected. ycMS 
is the beam rapidity in the collision center-of-mass frame. The two numbers of tracks/event refer to the strict and loose track 
quality cuts (see text). 



Beam energy (GeV) 


V^7(GeV) 


Vcms 


Year taken 


# of events 


Tracks/event 


20A 


6.3 


1.88 


2002 


140k 


58 - 66 


30A 


7.6 


2.08 


2002 


170k 


100 - 115 


40A 


8.7 


2.22 


1999 


160k 


141 - 164 


80A 


12.3 


2.56 


2000 


140k 


288 - 347 


158A 


17.3 


2.91 


1996 


120k 


455 - 568 



in a downstream calorimeter (VCAL) by projectile spec- 
tator nucleons. The geometrical acceptance of the VCAL 
calorimeter was adjusted for each energy using a collima- 
tor (COLL) pi]. 

Details of the data sets used in this analysis arc given in 
Table 1. We only included the 3.5% most central Pb+Pb 
collisions at each energy. The tight centrality selection 
was used to minimize event- by- event fluctuations due to 
the residual dependence of the kaon to pion ratio on sys- 
tem size and hence centrality. For the selected centrality 
range, we find that the remaining relative variation of 
the K/tt and (p + p)/n ratios with E veto agrees within 
1% for the different collision energies. In principle, this 
systematic variation could be washed out to varying de- 
grees, depending on the resolution of the VCAL at the 
various beam energies. To constrain this effect, we used a 
Glauber calculation in combination with a parametriza- 
tion of the E veto resolution of the VCAL calorimeter, 
taken from 



15] 



With these simulations, we found the 
RMS variation in the number of participating nucleons 
in the selected collisions at the lowest energy to be 3.9%, 
compared to 3.6% at the highest energy. As the abso- 
lute contribution of the K/tt centrality dependence in 
the 3.5% E veto bin on the final dynamical fluctuation 
&dyn is less than 1%, the additional energy dependence 
induced by the varying E veto resolution is expected to be 



less than 0.1%. No additional corrections are made for 
these contributions. 



B. dE/dx particle identification 

The overall sensitivity of the K/tt fluctuation measure- 
ments depends crucially on the resolution and stability 
of the event-by-event dE/dx particle identification. To 
optimize the stability of the dE/dx measurement with re- 
spect to time, variations of event multiplicity and possi- 
ble background contributions, only the energy loss of the 
track in the MTPCs was used in this analysis. The track- 
ing information from VTPC's was used in the momen- 
tum determination and for rejecting background from 
secondary interactions and weak decays. 

To eliminate a significant multiplicity dependence of 
the dE/dx due to the high charge load on the TPC read- 
out chambers in central Pb+Pb events, a correction algo- 
rithm developed specifically for the MTPCs was applied 
fltij ]. The response function of the TPC amplifier/shapcr 
was determined by running an iterative shape fitting al- 
gorithm over the MTPC data. Using the response func- 
tion obtained from this procedure, a channel-by-channcl 
correction of the raw TPC charge measurements is per- 
formed. This unfolding procedure takes into account the 
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charge history of sets of neighboring channels, which are 
coupled via crosstalk effects through the sense wires of 
each TPC readout chamber. 

These corrections improve the average dE/dx resolu- 
tion by about 30% from a^E/dx/ (dE/dx) = 5-6% to the 
final value of 3.9% for central Pb+Pb collisions. More 
importantly, the observed multiplicity dependence of the 
dE/dx measurement was reduced by more than 90%, 
leaving a change of less than 0.3% over the selected mul- 
tiplicity range. The values are quoted for the highest 
beam energy, corresponding to the highest multiplicities 
and charge loads in the TPCs. The performance of the 
algorithm was also cross-checked at each beam energy. 

The MTPCs predominantly detect particles in the for- 
ward hemisphere of the collisions. This requires particle 
identification to be performed in the relativistic rise re- 
gion, leading to a rather small separation in energy loss 
between different particle species. To allow for separa- 
tion of kaons and protons, a cut in total momentum of 
Ptot > 3 GcV/c is applied for all data sets. Above this 
cut, the achieved dE/dx resolution translates into an av- 
erage separation of pions from kaons of about 2.1 a and 
of kaons from protons of about 1.8 a. 

In combination with the geometrical acceptance of the 
TPCs and the track quality cuts, the p to t cut determines 
the acceptance function for particles entering the anal- 
ysis at each beam energy. To minimize the rapidity ac- 
ceptance shift for kaons, relative to CMS rapidity, the 
magnetic field settings for the two magnets were changed 
for each beam energy. Note that in a fixed target setup 
it is not possible to simultaneously keep the relative ra- 
pidity acceptance constant for multiple particle species 
with different masses. Figure [2] shows the phase space 
distributions of accepted particles at the different colli- 
sion energies, assuming pion, kaon and proton mass for 
the left, center and right columns. For each plot, the 
solid line indicates the default p tot > 3 GcV/c. The fig- 
ure illustrates the shifting acceptance, in particular for 
pions, in (pt, y) due to the changing acceptance and the 
constant PID cutoff in lab total momentum, ptot- Shown 
also are dashed lines for alternative ptot cuts that were 
applied to test effects of acceptance variations. 



III. ANALYSIS METHOD 
A. Event-by-event ratio determination 

Due to the finite number of particles detected for sin- 
gle Pb+Pb collisions and the limited particle identifica- 
tion capabilities of a TPC in the relativistic rise region, 
it becomes crucial to make maximum use of the avail- 
able information when constructing an estimator for the 
event-by-event particle ratios. The limited separation be- 
tween different particle species in the dE/dx measure- 
ment does not permit a simple counting of particles in 
this experiment. Instead, we use the unbinned distribu- 
tion of particles in (p, dE/dx) for each event to extract 
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FIG. 2: Phase-space distributions of charged particles used 
for the fluctuation analysis at all collision energies, from 
^/s NN = 6.3 to 17.3 GeV. Data are shown using pion, kaon 
and proton mass. The solid line indicates the 3 GeV/c cutoff 
in total lab momentum. The dashed line indicates cutoffs at 
each energy that were used for the study of acceptance effects 
(see text). 

just two parameters, K/ir and (p + p)/ir by p erforming 
an event-by-event maximum likelihood fit O El. 

At each energy, the probability density functions 
(PDFs) for the particle momenta, normalized to unity, 
F m (p), are determined for each particle species (to = 
kaons, pions, protons, electrons). We also evaluate the 
normalized PDFs for the truncated mean energy loss, 
f m (p', dE/dx), as a function of particle momentum for 
each species. 

For the event-by-event fit, the relative yield of different 
particle species is characterized by parameters m , such 
that O m = 1- These parameters arc determined for 
each event by maximizing the following likelihood func- 
tion: 

n 

L = []E ®mF m (pi)f m (pi; (dE/dx),)}. (1) 

i=l m 

Here pi denotes the observed momentum vector and 
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FIG. 3: (Top) Event- averaged dE/dx distributions (symbols) compared to the extracted dE/dx probability density functions 
(solid lines) for protons, kaons, pions and the sum of all particles. For all energies, data distributions and projected PDFs are 
shown for a narrow momentum range (Ap tot « 1 GeV/c) for 200 MeV/c < pr < 400 MeV/c, close to y = assuming kaon 
mass. (Bottom) Fit residuals for the sum of dE/dx probability density functions shown in the upper panels. 



(dE/dx)i the specific energy loss for each particle i in 
the event. The parameters G m were constrained to be 
positive for the maximum likelihood fit. The maximum 
likelihood fit then yields directly the event-by-event K/tt 
and (p + p)/n particle ratios. No other parameters are 
fitted event-by-event. 

The PDF used in the fit are determined from event- 
averaged data for each beam energy. The particles from 
each sample are split into bins in total momentum ptot, 
transverse momentum px, azimuthal angle <f> and charge. 
The dE/dx distribution of the particles in each bin is fit- 
ted by four Gaussian distributions, one for each particle 
species. The fit results for all bins are stored, providing 
the description of the PDF for each bin. The position and 
width of the Gaussian distributions describe the dE/dx 
PDF, f m , while the integrals of the distributions in the 
momentum bins describe the PDF of the momentum dis- 
tributions, F m . 

The top row in Figure [3] shows the event-averaged 
dE/dx distributions for positive particles in narrow mo- 
mentum slices Apt t « 1 GeV/c for each of the 5 collision 
energies. For each energy, the selected region was cho- 
sen close to rapidity y = 0, assuming kaon mass, with 
transverse momenta 200 MeV/c < pr < 400 MeV/c. 
Also shown in the top row as solid lines are the extracted 
PDFs for p, K, 7r + and e + in the selected momentum 
slice and the PDF sum. The bottom plots show the 
residual difference between data and sum PDF, illustrat- 
ing the quality of the dE/dx parametrization. A small, 
systematic structure in the residuals can be observed at 
all energies, resulting from small non-Gaussian tails in 
the dE/dx distributions. The dominant contribution to 



these tails comes from a variation in the number of hits 
in individual particle tracks, where tracks with a smaller 
number of hits lead to a wider dE/dx distribution. We 
did not include these non-Gaussian tails in the dE/dx 
PDFs. The resulting effect on the extracted fluctuations 
was tested by changing the dE/dx scaling of the width 
of the dE/dx fits (see below), which changes the relative 
contribution of the non-Gaussian tails, and by varying 
track quality cuts, and is included in the final systematic 
uncertainties. 

The results of the event-by-event fits for the five data 
sets are presented in Figure[31 showing the cvent-by-event 
ratio distributions for A"/7r (left) and {p + P)/tt (right). 
Points show the K/tt and (p+p) / 7r distributions for data, 
while the histogram shows the corresponding reference 
distributions obtained for mixed events (see below). 



B. Measurement of non-statistical fluctuations 

The next step in the analysis is to extract the strength 
of non-statistical event-by-event fluctuations from the ob- 
served distributions of A/7r and (j> + p)/ir. The relative 
width cr, defined as a = RMS/Mean x 100 [%], of the 
measured event- by- event particle ratio distributions can 
be decomposed into three contributions: 

1. due to the finite number of particles produced and 
observed per event, the ratio of particle multiplici- 
ties measured event- by- event will exhibit statistical 
fluctuations with a width dictated by the individual 
particle multiplicities within the acceptance, 
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events are constructed by randomly selecting measured 
particles from different events and combining them into 
artificial events, while reproducing the multiplicity dis- 
tribution of the real events. By construction, mixed 
events have on average the same particle ratios as the 
real events, but no internal correlations. In creating the 
mixed events, the dE/dx measurement for each particle 
is carried over to the mixed events. This allows us to esti- 
mate not only the finite-number fluctuations, but also the 
effect of the finite dE/dx resolution using mixed events. 
The mixed events are subjected to the same fit procedure 
as the real events and the RMS- width <J m i X of the mixed 
event particle ratio distribution is obtained. The results 
of the event-by-event fit on mixed events are compared 
in Figure Q] to those from real events. As the figure illus- 
trates, the contributions from the statistical and experi- 
mental fluctuations dominate the event-by-event particle 
ratio distributions, and have to be subtracted carefully. 

To extract the physically relevant non-statistical fluc- 
tuations from the data distributions as shown in Fig- 
ure |H we calculate for each energy the RMS widths of 
the event-by-event ratio distributions for data and mixed 
events. The non-statistical fluctuations Odyn are then es- 
timated by subtracting the relative RMS width of data 
and mixed event distributions in quadrature: 



<Jdyn = sign(a 2 dat 



mixed 



mixed I 



(2) 



It is important to note that non-statistical fluctuations 
7^ can arise from correlated particle production. 
Possible sources of particle correlations include energy- 
momentum and quantum-number conservation, quantum 
correlations and the decay of resonances (which includes 
strong and electromagnetic decays and a possible con- 
tamination due to weak decays). In particular, the defi- 
nition of a dyn 

allows for aj ata < o 2 mlxed , which could re- 
sult from correlated production of the two particle species 
forming the ratio, as expected for pion and proton pro- 
duction from hyperon decays. 



FIG. 4: (Left) Distributions of the event-by-event K/tt ratio 
for data (points) and mixed events (histogram), for all five 
collision energies. (Right) Distributions of the event-by-event 
(p+p)/n ratio for data (points) and mixed events (histogram), 
for all five collision energies. No acceptance corrections were 
applied to the particle ratios shown. 

2. due to non-ideal particle identification, these sta- 
tistical fluctuations will be smeared by the exper- 
imental dE/dx resolution and the event-by-event 
fitting procedure, 

3. due to genuine non-statistical fluctuations, which if 
they exist are superimposed on the background of 
statistical and experimental fluctuations. 

The combined contributions of finite number fluctuations 
in the particle multiplicities and effects of detector resolu- 
tion are estimated using a mixed event technique. Mixed 



C. Systematic uncertainty estimates 

As the comparison of the event-by-event particle ratio 
distributions for data and mixed events shows, the over- 
all width of the distributions is dominated by the fluc- 
tuations due to finite particle statistics and resolution. 
Clearly, the measurement relies on the absence of addi- 
tional detector-related fluctuations not accounted for by 
the mixed-event technique and on the reliability of the 
event-by-event fitting procedure. Several tests on modi- 
fied mixed events and on input from transport model cal- 
culations were performed to verify these requirements. In 
addition, we varied the event and track selection criteria, 
as well as the parameters of the dE/dx fits to estimate 
the associated systematic uncertainties. The results of 
these tests are described in the following section. Fig- 
ures [SHE] show the results for K/tt ratios. For all tests, 
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the uncertainties on (p + p)/ir fluctuations were found to 
be smaller or equal to those on K/tt fluctuations. This 
is expected, as protons are farther separated in dE/dx 
from the dominant pion contribution. 

Figure [9] summarizes the results of all of the system- 
atic tests described in this section. The error bars in that 
figure show the estimated statistical uncertainty on the 
change of the K/tx dynamical fluctuations for each varia- 
tion of fitting procedure or track and event selection. For 
tests performed on data, as discussed in Sections IIII CI 3 
and 6, MC estimates of the correlation coefficient for the 
fit results on different samples were used to take into ac- 
count the large overlap between the samples when e.g. 
removing the top or bottom 1% of events or changing 
track quality cuts. 



1. Linearity of event-by-event fit 

Using the PDFs as a function of (p, dE/dx) extracted 
from the event-averaged data, the particles in the mixed 
event pool can be assigned a particle species, such that 
on average they conform to the known PDFs. This as- 
signment then allows us to create modified mixed events 
with the desired multiplicity distribution and adjustable 
ratios of kaons, pions and protons. Thus we are able to 
vary both the mean and fluctuations in e.g. the K/n ra- 
tio in the modified mixed events. Using this technique, 
we created modified mixed events with K/n ratios vary- 
ing from 0.08 to 0.26 using the multiplicity distribution 
for y/s NN = 17.3 GeV and subjected these mixed events 
to the event-by-event fit. As shown in FigO this test 
demonstrated that over a wide range in K/n the event- 
by-event ratio estimator is linear with respect to the true 
K/n ratio in the event. A small bias of 0.002 in the ex- 
tracted particle ratio is observed. As the same fitting 
procedure is used on data and mixed events, this bias 
will not affect the estimate of <7d yn - 



2. Sensitivity of fluctuation measurement 

The modified mixed events discussed in the previous 
section also allow us to test the sensitivity of the ex- 
traction of non-statistical ratio fluctuations. Using the 
multiplicity distribution of the 17.3 GeV dataset, non- 
statistical Gaussian variations of magnitude <J™ yn in the 
K/n ratio were introduced as a "signal" in the creation 
of the modified mixed events, in addition to the statis- 
tical event-by-event variation. For values of u™ yn from 
to 20%, these modified mixed events were then run 
through the full analysis procedure, using unmodified 
mixed events as the statistical reference. The result of 
this study is shown in Figure [6l where the extracted fluc- 
tuations are compared to the known input amplitude. 
The extracted magnitude of fluctuations is found to be 
equal to the input amplitude within statistical errors, 



even for small input fluctuations of w 1%. A small sys- 
tematic uncertainty of 0.4% is assigned to the fitting pro- 
cedure for both Kj n and (p + p)/n non-statistical fluc- 
tuations. 
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FIG. 5: Difference between 
measured and input K/n 
ratio versus average input 
K/n ratio for mixed events, 
showing that the event-by- 
event fit is linear with re- 
spect to the input ratio. 



FIG. 6: Relative devia- 
tion of extracted K/n non- 
statistical fluctuations from 
the input value, a^ n , as 



a function of a 



dyn 



for 



modified mixed events (see 
text). The shaded band 
indicates the assigned esti- 
mate of systematic uncer- 
tainty from the fitting pro- 
cedure as a function of am- 
plitude. 



Another concern is the highly non-Gaussian nature 
of the event-by-event K /it distributions for low particle 
multiplicities. To check the validity of equation 2 in this 
case, the observed distributions for mixed events were 
fitted with Gamma-distributions at each energy. The 
Gamma-distributions were folded with Gaussians corre- 
sponding to the magnitude of non-statistical fluctuations 
at each energy and sampled to reproduce the input his- 
tograms for <Tj ata . The unfolded Gamma-distributions 



were sampled to obtain cr^ lfaed . The resulting a^ yn was 
plotted for each energy, in comparison to the correspond- 
ing input value, in Figure |9ja) . This study demonstrates 
that the magnitude of the Gaussian input fluctuations is 
recovered from applying equation 2 to the RMS width of 
the non-Gaussian signal and reference ratio distributions. 
The effect of low particle multiplicities was also tested 
by repeating the analysis for the highest energy dataset 
while randomly removing particles from each event, to 
resemble the multiplicity distributions in the lower en- 
ergy datasets. The extracted particle ratios were found 
to be independent of the fraction of removed particles 
within statistical errors. 

Finally, the analysis method was tested on events sim- 
ulated with the UrQMD transport code For these 
events, non-statistical fluctuations as a function of beam- 
energy were determined using two methods. The first 
method determined the event-by-event particle ratios by 
counting the different particle multiplicities in the NA49 
acceptance. The second method assigned a dE/dx value 
to each particle according to the (p, dE/dx) PDFs and 
then performed an event-by-event fit to extract the par- 
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ticle ratios. For each method, the non-statistical fluc- 
tuations were extracted by comparison to the respective 
mixed event reference. As Figure Mb) demonstrates, the 
two methods yield the same magnitude of non-statistical 
fluctuations at each collision energy, within statistical er- 
rors. To account for a possible small energy dependence, 
the systematic uncertainty shown by the shaded area in 
Figure [Ub) was assigned. 




0.05 0.1 °% 0.05 0.1 

Event-by-event KJ% Event-by-event K/ji 



3. Fluctuations in the tails of ratio distributions 

In Figure [H several outlier events in the ratio distri- 
butions can be seen, in particular for the higher collision 
energies. To test whether the observed non-statistical 
fluctuations are due to these outliers or fluctuations in 
the tails of the ratio distributions, we removed 1% of 
events with the highest K/tt ratio in data and mixed 
events at each energy. The resulting change of the ex- 
tracted fluctuation signal is shown by the full symbols 
in Figure Etc) . Also plotted in this figure is the result of 
removing the 1% of events with the lowest K/tt ratio. As 
expected, removing either of the tails of the distribution 
leads to a slight reduction in the extracted dynamical 
fluctuations ad yn - In particular for removing the high 
tail of the distribution, the change is small for most en- 
ergies. Removing the low end of the distribution has a 
10% effect at the low energies. Smaller changes were seen 
for a corresponding study of (p + p)/tt fluctuations. 

The K/tt ratio distributions for the lowest two ener- 
gies develop a spike at very small K/tt ratios. Figure [7| 
shows an expanded view of the small K/tt ratio region 
for the ^/s^T = 6.3 GeV data set. We examined the 
events in this region for signs of detector or reconstruc- 
tion failures, but found no anomalies compared to the 
bulk of events. The spike is also present in the mixed 
event sample, suggesting that it arises from a combina- 
tion of low kaon multiplicities and the finite resolution 
of the event-by-event maximum likelihood fit. It is im- 
portant to note that the fit does not allow negative K/tt 
values. To confirm this hypothesis, we studied the ra- 
tio of the event-by-event K/tt distributions in data and 
mixed events, shown in Figure [H] This ratio does not 
exhibit a spike at small K/tt ratios. The shape of this 
ratio is consistent with the presence of dynamical fluc- 
tuations in the data events and with the observed small 
decrease in (Tdyn when removing events in the low tails of 
the distributions from the analysis. 

Finally, we performed further simulations using the 
UrQMD generator. For generated UrQMD events, event- 
by-event K/tt distributions were obtained by counting 
particles in the NA49 acceptance and alternatively by 
assigning energy loss information to each particle and 
performing the same event-by-event likelihood fit as de- 
scribed before. As expected, a broadening of the event- 
by-event K/tt distribution is seen when using the like- 
lihood fit as compared to particle counting. This leads 
to an enhanced number of events in the K/tt ss bin, 



FIG. 7: Expanded view FIG. 8: Ratio of the event- 

of the low K/tt ratio re- by-event K/tt ratio distri- 

gion of the event-by-event bution for data and mixed 

K/tt ratio distribution of events of the ^/s NN = 

the y?s NN — 6.3 GeV data 6.3 GeV data set. 
set for data (points) and 
mixed events (histogram). 



reproducing what is seen in the data. However, this en- 
hancement is again well reproduced in the corresponding 
mixed events and no bias is found compared to an analy- 
sis of dynamical fluctuations based on direct counting of 
particles. 

In summary, our studies show that the extreme tails 
of the distribution do not dominate the overall fluctua- 
tion strength. In particular the shape of the K/tt distri- 
butions for small K/tt is well reproduced by the mixed 
event procedure and in UrQMD studies using simulated 
energy loss information. 



4- Uncertainty in dE/dx scaling 

In determining the dE/dx PDFs, we assume that the 
width of the dE/dx distributions for the different par- 
ticle species scales with the average dE/dx as dE/dx a , 
with a — 0.65 obtained in previous studies of dE/dx par- 
ticle identification [2C| . As a function of momentum p, 
the width for a given particle species also changes due to 
variations in the observed TPC track length and due to 
variations in the local track density. In combination with 
the small separation between different particle species, 
this makes a precise determination of a difficult. To test 
the effect of the uncertainty in a on the observed fluctu- 
ations, we repeated the analysis, starting from the deter- 
mination of PDFs, for a much larger value of a = 1.15. 
The difference between the fluctuation results obtained 
for a = 1.15 and a = 0.65 is shown in Figure [9^d) . The 
fluctuations are seen to be rather robust against changes 
in a, with the difference only reaching up to 1% for the 
lower collision energies. Even smaller variations are seen 
in the (p + p)/tt fluctuations. The actual uncertainty in 
a is estimated to be only ±0.15, leading to the system- 
atic uncertainty estimate shown as the shaded band in 
Figure [^d). 
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5. Track quality cuts 

In order to test for possible systematic distortions of 
the fluctuation measurement not included in the MC sim- 
ulations, we used two sets of quality cuts in the event 
and track selection for the analysis. The correspond- 
ing datasets should have different sensitivity to detector 
and tracking effects, such as contributions from weak de- 
cays and secondary interactions, or effects from varying 
track densities on reconstruction and particle identifica- 
tion. The quality cuts are based on the distance of closest 
approach of the extrapolated particle trajectory to the 
main vertex, the length over which the track was mea- 
sured in the TPC and the number of measured points on 
the track. For the set of loose quality cuts, the cuts on the 
extrapolation of tracks to the primary vertex and a cut on 
the track fit quality were removed. As shown in Table 1, 
the difference in multiplicity for the two samples varies 
with increasing track density and changing acceptance, 
as the collision energy increases, and ranges from 10% 
at the lowest energy to about 20% at the highest energy. 
The track selection using loose track quality cuts yields 
more tracks per event to include in the event-by-event 
particle ratio estimation, but possibly also increases the 
contamination of the track sample with non-primary and 
background particles. 

The final value of the fluctuation signal presented be- 
low is calculated as the arithmetic mean of the results 
from both samples. Figure[9Je) shows the deviation from 
the average result for the tight track-cut sample. The 
shaded area indicates the estimated systematic uncer- 
tainty in the final result for K/tt fluctuations related to 
the track selection. The same study was performed for 
(p+p)/ir fluctuations, resulting in a smaller variation 
than seen for K/tt. The corresponding contribution to 
the systematic uncertainty on the (p+p) / tt non-statistical 
fluctuations is shown by the dashed lines in Figure [He). 



6. Effects of acceptance variation 
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FIG. 9: Results of various studies on systematic uncertainties 
in the determination of o-d yn (K/n): (a) extraction of adyn 
using equation 2. (b) dE/dx fit versus particle counting in 
UrQMD. (c) removal of tails of event-by-event distributions, 
(d) variation of dE/dx width scaling, (e) variation of track 
quality cuts, (f) variation of ptot cutoff. Shaded bands indi- 
cate the contribution to the systematic error for K/tt fluctua- 
tions. Where applicable, dashed lines indicate corresponding 
contribution to systematic error for (p + p)/n fluctuations. 
See text for details. 



As discussed above, the constant PID cutoff in lab to- 
tal momentum p tot > 3 GeV/c, below which particles 
species cannot reliably be separated, leads to a signif- 
icantly changing acceptance at different beam energies 
(see Fig. 2). To study the influence of the changing ac- 
ceptance relative to the lowest beam energy, we modified 
the cutoff for the four higher beam energies to 3.6, 4.4, 
6.3 and 9.1 GcV/c, respectively These higher cutoffs re- 
move low pt particles near midrapidity approximating 
the default cutoff for the ^/s NN = 6.3 GeV data sam- 
ple. The resulting change in the K/tt fluctuation sig- 
nal is plotted as the filled symbols in Figure Eff), show- 
ing a rather small difference compared to the standard 
acceptance cut. We also studied a reduced p tot cut of 
1.45 GeV/c for the lowest energy data set. The region 
1.45 GeV/c < ptot < 3 GeV/c yields little separation 



between kaons and protons, but extends the pion accep- 
tance closer to mid-rapidity. The corresponding fluctua- 
tion signal is shown as the open symbol in Figure [^f). 
Again, no strong change in the K/tt non-statistical fluc- 
tuations is observed. Note that the event multiplicity is 
changed by up to 30% by the modified cuts, such that 
part of the change observed is likely due to statistical 
fluctuations. As any remaining variation could be due 
to the physical nature of the non-statistical fluctuations, 
no systematic uncertainty is assigned for the acceptance 
variations. 

While this study indicated that the results for K/tt 
fluctuations only have a moderate dependence on the de- 
tails of the acceptance, we need to stress again that a de- 
tailed comparison to theoretical models needs to take the 
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experimental acceptance into account. As the measure- 
ment is performed within a limited acceptance, particle 
ratio fluctuations within the acceptance can in principle 
be caused by non-statistical variations of e.g. the kaon 
phase space distribution, rather than a non-statistical 
variation of the K/tt ratio in full phase space. Note how- 
ever that purely statistical variations of the phase-space 
distributions are removed by comparison to the mixed 
event reference. 



IV. RESULTS AND DISCUSSION 

The event- by-event K/tt and (p + p)/tt ratios shown in 
Figure [4] exhibit smooth distributions at all five collision 
energies and show no evidence for the existence of distinct 
event classes. Qualitatively the K/tt ratio distributions 
are found to be wider than the statistical reference for 
all energies, while the (p + p)/n ratio distributions are 
narrower than the respective reference. 



The deviation from the statistical reference 



as 



defined in equation^ is plotted in FigurefTOlfor K/tt fluc- 
tuations (top) and (p + p)/tt fluctuations (bottom) . As 
stated in Section IIII C 5| the final data represent the av- 
erage of the loose and tight track quality cut results. The 
observed non-statistical fluctuations of the K/tt ratio are 
positive and decrease with collision energy. For the high- 
est collision energy, a value of <Jd V n — 3.2 ± 0.4(stat.) ± 
0.5(syst.)% is observed, in agreement with the results re- 
ported in [T3| for an independent dataset. Towards lower 
collision energies, a steep increase of the fluctuation sig- 
nal is observed, the K/tt fluctuations reaching a value of 
Pdyn = 7.9±1.2(stat.)±1.0(syst.)% for the lowest energy. 
The increase is qualitatively consistent with calculations 
assuming the presence of a deconfinement phase transi- 
tion at the lowest collision energies of the CERN SPS [2lj . 
This scenario predicts small and energy-independent fluc- 
tuations in the K/tt ratio in the QGP phase, while for 
the confined phase large fluctuations are expected that 
increase towards lower collision energies. 

For the {p + p)/tt ratio, which has not been previously 
reported, the non-statistical fluctuation measure <Jd yn is 
negative for all collision energies. For the lowest collision 
energy a value of (Jdyn = — 8.1 ± 0.4(stat.) ± 1.0(syst.)% 
is observed. The negative value of (Jdyn for (p + P) /tt 
fluctuations reflects the fact that the width of the data 
distribution is smaller than the width of the distribution 
of mixed events for this ratio. A negative fluctuation sig- 
nal of the event-by-event (p + p) /tt ratio can result from 
resonance decays into pions and protons. Such contribu- 
tions correlate the pion and proton multiplicities event- 
by-event, and thus lead to smaller fluctuations in the ra- 
tio than expected from independent particle production. 
The magnitude of the negative fluctuation signal in the 
(p + P)/tt channel is then related to the relative contri- 
bution of resonance decay products in the final state of 
the collision. 

In order to further estimate the influence of known cor- 
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FIG. 10: (Color online) Energy dependence of the event-by- 
event non-statistical fluctuations of the K/tt ratio (top panel) 
and the (p + p)/7r ratio (bottom panel). Filled symbols show 
data, open symbols show calculations with the UrQMD trans- 
port code, using NA49 acceptance tables. Systematic uncer- 
tainties are shown as brackets. Data of the STAR collabora- 
tion are taken from 12411 . 



relations in the events, like those induced by conservation 
laws and particle production from resonance decays, we 
analyzed fluctuations of the K/tt and (p + p) /tt ratios in 
a string-hadronic cascade model, UrQMD [lj|. In this 
model, by construction, no fluctuations due to a poten- 
tial phase transition are present. For this study, large 
samples of UrQMD events were generated at all five col- 
lision energies and then subjected to an acceptance filter 
modeling the NA49 detector system. The accepted final 
state particles were counted and the corresponding ratios 
were formed for each event. The energy dependence of 
the event-by-event (p+p)/tt ratio fluctuations in UrQMD 
closely matches the energy dependence observed in the 
data, as shown in the bottom plot in Figure [TO] This 
lends further support to interpreting the negative fluctu- 
ation signal as resulting from resonance decays, which are 
the dominant source of {p+P)/tt fluctuations in UrQMD. 

The energy dependence of the K/tt ratio fluctuation 
signal however is not reproduced in the cascade model, 
which gives an energy independent fluctuation signal. 
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Further studies are needed to determine whether the fi- 
nite non-statistical K/tt fluctuations in UrQMD originate 
from resonance decays or correlated particle production 
due to conservation laws. 

The K/tt fluctuation signal in the data taken at 
y/s NN = 17.3 GeV collision energy was found to be 
consistent with calculations performed assuming a grand 
canonical ensemble without enforcing local conservation 
laws (22)]. The small fluctuations at this energy are thus 
consistent with every event being a random sample from 
the identical thermal ensemble. The interpretation of the 
observed increase of K/tt fluctuations towards lower en- 
ergies is complicated by the simultaneous decrease in true 
particle multiplicities towards lower energies. Stephanov 
has pointed out that the observed energy dependence 
could approximately be understood as the combination 
of a fixed pair-wise correlation strength, in combination 
with the known dependence of the average event multi- 
plicity on collision energy (23[. Within the experimental 
resolution of the present measurement alone, the multi- 
plicity scaling expected from this argument can not be 
ruled out. 

Recently, the STAR collaboration has reported pre- 
liminary results for non-statistical K/tt fluctuations in 
central Au+Au collisions in the RHIC energy range from 
v /s^7 = 19.6 to 200 GeV [Hj]. The reported results for 
19.6 GeV are in agreement with the values quoted here 
for y/s NN = 17.3 GeV, within statistical errors. Although 
a detailed quantitative comparison of the results suffers 
from the differences in the NA49 and STAR acceptances, 
it is important to note that the STAR data show no sig- 
nificant energy dependence, in contrast to the results re- 
ported here for the SPS energy range. If confirmed, this 
would rule out, or at least limit the range of validity, for 



the proposed multiplicity scaling discussed above. 



V. SUMMARY 

We have presented a study of the energy dependence 
of non-statistical event-by-event fluctuations of the K/tt 
and (p + p)/tt ratios in the energy range 6.3 < y/s NN < 
17.3 GeV. A strong increase of the K/tt fluctuation signal 
is observed towards the low end of this energy range. The 
increase of the signal is not reproduced by the UrQMD 
hadronic cascade model, suggesting the onset of a new 
source of fluctuations. The same model is however able to 
describe the strong fluctuations seen in the (p + p) /tt ra- 
tio. The domain of increased K/tt fluctuations coincides 
with the energy range where a change in the behavior 
of event-averaged hadronic observables in heavy-ion col- 
lisions occurs, relative to elementary collisions. Further 
theoretical calculations will be needed to evaluate the rel- 
evance of this observation for a possible interpretation in 
the context of a deconfinement phase transition. 
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